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ABSTRACT 

Context. Both remote-sensing measurements using the interplanetary scintillation (IPS) technique and in situ measurements by the 
Ulysses spacecraft show a bimodal structure for the solar wind at solar minimum conditions. At present what makes the fast wind fast 
and the slow wind slow still remains to be answered. While a robust empirical correlation exists between the coronal expansion rate 
f c of the flow tubes and the speeds v measured in situ, further data analysis suggests that v depends on more than just f c . 
Aims. We examine whether the non-radial shape of field lines, which naturally accompanies any non-radial expansion, could be an 
additional geometrical factor. 

Methods. We solved the transport equations incorporating the heating due to turbulent Alfven waves for an electron-proton solar wind 
along curved field lines given by an analytical magnetic field model, representative of a solar minimum corona. 
Results. The field line shape is found to influence substantially the solar wind parameters, reducing the asymptotic speed by up to 
~ 130 km s -1 or by ~ 28% in relative terms, compared with the case neglecting the field line curvature. This effect was inteipreted in 
the general framework of energy addition in the solar wind: Relative to the straight case, the field line curvature enhances the effective 
energy deposition to the subsonic flow, resulting in a higher proton flux and a lower terminal proton speed. 

Conclusions. Our computations suggest that the field line curvature could be a geometrical factor which, in addition to the tube 
expansion, substantially influences the solar wind speed. Furthermore, at solar minima although the field line curvature unlikely 
affects the polar fast solar wind, it does help make the wind at low latitudes slow, thereby helping better reproduce the Ulysses 
measurements. 
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1. INTRODUCTION 

Around solar minima the global distribution of the solar wind 
exhibits a simple bimodal structure: a band around the helioe- 
quator with an angular width ~ 50° of the slow wind (speed 
v < 500 km s ) is sandwiched betwee n the high-latitude fas t 
winds (v > 650 km s" 1 ) (e.g., Plate 1 in lMcComas et al.ll2000l) . 
This picture was established indirectly by interplanetary scin- 
tillation (IPS) measurements in 1974 and 1985, correspond- 
ing to the minimum of solar cyc les (SC) 20 and 21, respec- 
tively dKoiima & Kakinumalll99Ch . However, it became defini- 
tive only when Ulysses finished th e first polar orbit whic h en- 
compasses the minimum of SC 22 dMcComas et al.ll2000l) . The 
bimodal picture continued for the SC 23 minimum as shown 
by Ulysses durin g its third polar orbit dMcComas et al.l 12.0081 
lEbert et al.ll2009l 

What makes the fast solar wind fast and slow wind slow? 
In the context of Ulysses measurements, what produces the 
pronounced latitudinal variation of v? In light of the well- 
established empirical relat ion between y and the c oronal expan- 
sion rate f c of flow tubes dWang & Sheelevlll990l) . one may ex- 
pect that the speed profile results from the differences in f c of 
tubes reaching different latitudes in interplanetary space. Models 
that incorporate heating rates closely related to the magnetic 
field B indeed suggest that varying the spatial distribution of B, 
or equivalently that of the expansion rate, can result in consid- 
erable changes i n wind paramet ers (for a summary of geomet- 
rical models, see ICranme r 2005; f or more rec ent models adopt- 
ing empirical heating terms, see IWangll20TTl) . For instance, the 



most recent models based on the anisotropic turbulence, which 
accommodate self-consistently the heating due to interactions 
of counter-propagating Alfven waves, indicate that by varying 
the B profile alone, the modeled latitudinal dependence of wind 
speed at 1 AU c an agree rather well with the Ulysses measure- 
ments (Fig. 12 in lCranmer et al . 2007|). 

While the f c - v relationship is strikingly robust, the signifi- 
cant spread in y in a given f c interval means that v is a function 
of more than f c , and the angular distance 9b of the field line foot- 
point to its nearest coronal hole boundary was pr oved a good 
candidate of this additional factor ( lArge et al.ll2004l) . For a given 
f c , the smaller 9b is, the lower the speed. As one naturally expects 
that the open field line becomes increasingly curved while its 
footpoint becomes closer to the coronal hole boundary, a natural 
alternative to 9b would be the shape of field lines. Surprisingly, 
the possible effect of field line curvature has not been examined. 
Instead, the available turbulent-heating-based models tend to as- 
sume a radially directed flow tube, even though the non-radial 
tube expansion has been examined in considerable detail (e.g., 
IChen & Hdl2002t ICranmer et al.ll2007HWangil201 ll) . 

The aim of this study is to examine how the field line shape 
affects the solar wind parameters in general, and whether it helps 
reproduce the Ulysses measurements in particular. It should 
be noted that while fully multi-dimensional models incorpo- 
rating turbulent or e mpiri c al heating rates ( e.g.. lUsmanov et aD 
' J00f IChen & Hul 120011: iLi et alJ [2004 ICohen et al.l J2007h 
Lionello et alj|2009t iFeng et alj|2010t Ivan der Hoist etallboiol) 



certainly have taken account of both a non-radial expansion and 
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the field line curvature, the two effects cannot be separated. 
Therefore we choose to solve the transport equations along field 
lines in a prescribed magnetic field configuration typical of the 
solar minimum corona, thereby allowing us to switch on and off 
the effect of field line shape. The solar wind model we employ 
incorporates a physically motivated boundary condition, distin- 
guishes between electrons and protons, adopts reasonably com- 
plete energy equations including radiative losses, electron heat 
conduction, and the solar wind heating from parallel propagat- 
ing turbulent Alfven waves. 

The paper is organized as follows. In section [2] we give a 
brief overview of the model and describe how it is implemented. 
The numerical solutions are then presented in section [3] Finally, 
section |4] summarizes the results, ending with some concluding 
remarks. 



2. Model Description and Solution Method 

We consider a time-independent electron-proton solar wind 
flowing in the meridional plane and incorporating turbulent heat- 
ing from Alfven waves. Since the equations are quite well- 
known, we simply note that the mass or energy exchange be- 
tween neighboring flow tubes is prohibited, permitting the de- 
composition of the vector equations (e.g., Eqs.(2) to (7) in 
iLi et al.ir2005l and references therein) into a force balance con- 
dition transverse to B and a set of conservation equatio ns along 
it. For mathematical details, please see the Appendix of ILi & LJ 
d2006l) . In this study we replace the former with a prescribed 
magnetic field, and treat in detail only the conservation equa- 
tions. The Alfven waves are assumed to propagate along B, and 
to be dissipated at the Kolmogorov rate gwav = poV /L c where 
p = nm p is the mass density with n and m p being the number 
density and proton mass, respectively. Moreover, Sv is the rms 
value of the wave velocity fluctuations, and L c the correlation 
length. By standard treatment, £? wav is entirely deposited to pro- 
tons, and L c - L c q(Bo/B) 1 ^ 2 where the subscript denotes val- 
ues at the field line footpoint at the Sun. This form of dissipation 
and its variants, first proposed by iHollwed (119861) . have been a 
common practice in solar wind mode ling ever since (see e.g., the 
review bv lHollweg & Isenberg|l2002l) . 

The conservation equations along an arbitrary magnetic line 
of force read 



(nvA)' 



= 0, 



p p \ r ) p 



P P 

r; + ( r -i)r,[in(vA)]' 
r -i 



nkgv 
(F W A)' 



6 -^ + Q s -L- { ML 
5t ^ A 



+ vF 



(i) 

(2) 

(3) 
(4) 



where v is the wind speed, A oc l/B is the cross-sectional area 
of the flux tube, and the prime ' denotes the derivative with re- 
spect to the arc length /, measured from the footpoint of the field 
line. The species pressures p s (s = e, p) are related to species 
temperatures T s via p s = nksTs in which kg is the Boltzmann 
constant. The gravitational constant is denoted by G, and M 
the mass of the Sun. Note that the gravitational acceleration is 
put in the form of a directional derivative. In addition, the force 
density acting on the fluid is due to Alfven waves, F = -p' w 
where p w = pdv 2 /2 is the wave pressure. The adiabatic index 



y — 5/3, and 6E s /5t represents the energy exchange rates of 
species s due to Coulomb collisions with the other species. To 
evaluate them, we take the Coulomb logarithm to be 23. For sim- 
plicity, the proton heat flux is neglected (q p = 0), whereas the 
electron one is taken to be the Spitzer law q e = -KT^ 2 T' e where 
k = 7.8 x 10~ 7 erg K~ 7/2 crrT 1 s _1 . Moreover, the heating rates 
applied to species s is denoted by Q s . For the protons, one has 
Qp = 2wav On the other hand, optically thin radiative losses are 
accounted for by the term L which is included in the electron ver- 
si on of Eq.Q, and des cribed by the parametrization scheme due 
to iRosner et al.l d 19781) . Some electron heating is also included, 
expressed by Q e = Q e exp(-l/ld) where Q e is the magnitude 
and Id the scale length of this heating. For simplicity, here we 
assume that this electron heating is associated with no momen- 
tum deposition. Equation describes the wave evolution, in 
which F w = p w (3v + 2va) is the wave energy flux density with 
va — Bj ^4-np the Alfven speed. 

Given the governing equations, one may say a few words on 
the effect of the magnetic field line shape before actually solving 
them. Consider the case, referred to as Case S, where we put dl = 
dr in ' (note that ' is also involved in q e ) and I = r-R G in Q e . Here 
r is the heliocentric distance and R Q the solar radius. As such, 
Case S differs from the original case, called Case C, only in that 
the field lines are straight in the former, but allowed to be curved 
in the latter. It then follows that if the wind energetics is not self- 
consistently considered, say, if one prescribes an r-distribution 
of the temperatures T e and T p , then Eqs.(Q]i and (f2]i indicate that 
no matter how curved the field line may be, the same r-profiles 
for n and v result. This applies even if the waves are considered 
as long as the wave dissipation is neglected. To understand this, 
one may r eadily derive f rom Eq.(H) a conservation law which 
reads (e.g. lHu et alll997l) 



(ojSA)' 
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where aiS = 2p H ,v^(l + Ma) is the wave action flux density, and 
Ma = v/\'a is the Alfvenic Mach number. If gwav = then one 
finds that uSA is a constant, enabling the wave pressure p w to 
be expressed as a functional of flow parameters as well as B. As 
such multiplying Eqs.dTJ and © by a factor dl/dr does not in- 
validate the equal signs. This means any two field lines that have 
identical radial distributions of B, T e and T p will yield identical 
radial distributions of n, v and p w , provided that the base val- 
ues for n and 6v are the same for the two lines. However, if the 
energy equation is explicitly considered, as we will do, and the 
adopted heating rates play a role in determining the solar wind 
mass flux, then the field line shape will likely play a role (see 
also the discussion section). 

For the merid i onal magnetic field, we adopt the 
iBanaszkiewicz et ail dl998h model for which we use the 
reference values [Q,K,ai] = [1.5, 1, 1.538]. This set of param- 
eters is supported by images obtained with the instruments on 
board the SOlar and Heliospheric Observato ry (SOHO) space- 
craft at solar minimum (see e.g., Fig.l in iForsvth & Marschl 
1 19991) . a nd has found extensive applications in solar wind mod- 
els (e .g., lHackenberg et al.ll2000l: ICranmer et al"1l2007t ILi & Lil 
120081) . Another parameter M which controls the magnetic 
field strength is chosen such that B at the Earth orbit is 3.5y , 
con sistent with U lysses measurements (ISmith & Baloghl[l995l: 
also lSmithll20TTl) . 

The meridional magnetic field configuration is depicted in 
FigHk where the thick contours indicate the lines of force along 
which we compute the expansion factor / = (A/Ao)(R /r) 2 , 
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Fig. 1. Adopted meridional magnetic field configuration in the corona, (a) A quadrant is shown where the magnetic axis points 
upward, and the thick contours delineate the lines of force for which the expansion factors / = (A/Ao)(R Q /r) 2 are plotted in (b) as 
a function of heliocentric distance r. The field lines are labeled with the colatitude at which they intersect the Earth orbit (1 AU). 
Here A denotes the cross-sectional area of the flow tube, the subscript represents values at the inner boundary R Q . 



given in FigQJ). The field lines are labeled with the colatitude 
at which the lines are located at 1 AU. From FigOJ one can 
see that the open field lines occupy only a minor fraction of 
the solar surface, with colatitudes 9 < 30° to be spec ific. If 
identi fying this area as the coronal hole (CH) (see e.g., I Wang! 
l2oTTh ." one can see that indeed when their footpoints approach 
the CH boundary, the open field lines become progressively 
more curved. Furthermore, Fig.QJ shows that while in the inner 
corona, say at r = 1.7 R , the expansion factor / increases when 
the field lines become closer to the CH edge, the opposite is true 
at larg e distances, say r > 20 R Q . Actually as was first pointed 
out bv lWang & Sheelevld 19901) . this tendency is due to the inclu- 
sion of the current sheet component in the magnetic field, and 
the low densities in the fast wind in interplanetary space were 
attributed to the the fact that flow tubes originating near the CH 
center underwent the most significant net expansion. 

Given a B profile along a field line, the governing equa- 
tions are put in a time-dependent form, discretized onto a non- 
uniformly spaced grid, and then advanced in ti me from an ar- 
bitrary initial state by a fully implicit scheme dHu et al. I fl997h 
until a steady state is reached, which is then taken as our so- 
lution. For the solutions found, the mass, momentum, and en- 
ergy are found to be conserved to better than 1%. At the top 
boundary (1 AU), all the flow parameters are linearly extrapo- 
lated for simplicity. At the base (1 R Q ) the flow speed v is de- 
termined by mass conservation, whereas the species tempera- 
tures are fixed at T e o = = To = 5 x 10 5 K, corresponding 
to mid-Transition Region. Because the processes taking place 
below this level should influence the mass supply to the sola r 
wind (e.g., lHansteen & Leerl Il995t iLie-Svendsen et al ] l2002h . 
they should be taken into account and can be well repr esented by 
the Radiation Equilibrium B oundary (REB) condition (Withbroe 
Il988t iLionelloet alJl200lh . The REB condition translates into 



«o = 4.8 x lQ 3 q e o in cgs units for our chosen Tq. As such 
the base density no is not fixed but allowed to adjust to the 
downward heat flux, which in turn depends on how the fluid 
is heated close to the base. Without the basal electron heating, 
the T e gradient and hence the electron heat flux density q e o at 
the base are too small, leading to HqTo <sz 1.6 x 10 14 cm" 3 K, 
the value inferred from measurements with the Solar Ultraviolet 
Measurements of Emitted Radiation ( SUMER) on SOHO for 
coronal holes dWarren & Hass ler 1999). The values adopted for 
Q e are Q e = 1.7 x 10" 5 erg cm" 3 s _1 , and l d = 0.06 R . With 
such a short scale length, Q e rapidly gives way to the more grad- 
ual proton heating Q wav for which we use a base wave amplitude 
dvo = 27 km s" 1 and a base value for the correlation length L c o of 
1.5 x 10 4 km. (Note that even though this L c o translates into ap- 
proximately 0.02 R Q , it does not mean that the wave dissipation 
is extremely concentrated at the base. A length scale comparable 
to Id involved in Q e can be given by I2wav/Gwavl> f° r which one 
finds a value ranging from 0.086 R Q at the base to 0.59 R at 
2 R for the solution labelled case C in Figj2j. 

3. NUMERICAL RESULTS 

To examine the extent to which the wind parameters are affected 
by the field line shape, we will compare solutions incorporating 
the curvature effect (Case C) and solutions neglecting it (Case 
S). 

Figure [2] presents the radial profiles of (a) the wind speed 
v and (b) the electron and proton temperatures T e and T p for 
the solutions computed along the tube labeled 88 in FigQJ. For 
such a field line the effect of the field line curvature is almost 
maximal. The solid and dashed curves correspond to Cases C 
and S, respectively. The asterisks in Fig|2^ refer to the sonic 
point where v equals the effective sound speed cj given by (e.g., 
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Fig. 2. Alfven wave-turbulence-based solar wind models incor- 
porating effects of field line shape. Radial profiles are shown 
for (a) the solar wind speed v, and (b) the electron and proton 
temperatures T e and T p . The curved (straight) case in which 
the field line shape is considered (neglected) is given by the 
solid (dashed) lines. In (a), the asterisks indicate the sonic point 
where v equals the effective sonic speed c$ defined by Eq.©, 
while the vertical bars represent wind speeds measured by track- 
ing small density in homogeneiti e s (blo bs) in SOHO/LASCO 
images as given by IWang et al.l d2000h . The triangles in (b) 
give the Hydrogen kinetic temperatures derived from H I Lya 
measure ments with SOHO/UV CS for a typical solar minimum 
streamer (IStrachan et al.ll2002h . 



Eq.(36) in lHackenberg et al.ll2000h 



kB_ 
m„ 



{Te + T p ) 



1 +3M A 



P 2(1 + Ma) 
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Moreover, in Figf2^ the vertical bars give the range of 
wind speeds derived by tracking a collection of small in- 
homogeneities (the blobs) in images obtained with the 
Large A ngle Spectrometric COronagraph (LASCO) on board 
SOHO dWang et al.l l2000h . The triangles in Figf2p repre- 
sent the Hydrogen kinetic temperatures derived from the H 
I Lya profiles measured with the Ultra Violet Coronagraph 
Spectrom eter (UVCS) on SO HO for a typical solar minimum 
streamer (IStrachan et al.ll2002h . 

The solution for Case C (solid curves) is in good agreement 
with available measurements. For instance, the model yields at 
1 AU the following parameters: Ve = 333 km s , (tiv)e = 3.65 x 
10 8 cm 2 s~ l ,T e = 1.4xl0 5 Kandr p = 2.7xl0 4 K, all of which 
match the in situ measurements by both H e lios (|Schwennll 19901) 
and Ulysses dMcComas et al.ll2000L 120081: lEbert et al.ll2009l) of 
the low latitude slow winds. Here the subscript E denotes the 
Earth orbit. Figure [2^ indicates that the wind speed starts with 



6.4 km s at the base, in line with th e Doppler speed mea- 
sured with Ne VIII (e.g. JXia et al.ll2003T) whose formation tem- 
perature is close to our base temperature. The speed profile is 
not monotonic but possesses a local minimum at 1 .63 R Q where 
v = 25.1 km s . The sonic point, denoted by rc, is located at 
4.15 R Q . Between 3 and 20 R Q the computed speeds are in line 
with the blob measurements. Moreover, for r < 3 R & the mod- 
eled v is no higher than 98.8 km s _I which w ould be below th e 
sensitivity of the Doppler dimming of H I Lya dKohl et alj2006h . 
and in this sense compatible with the fact that there is not much 
empirical knowledge about the proton speed in the inner portion 
of the slow wind. Moving to Figf2j5 one may see that due to the 
applied electron heating both T e and T p increase near the base. 
At 1.19 R e the electron temperature T e attains its maximum of 
9.8 x 10 5 K. Below this height, T e and T p stay close to each other 
due to Coulomb coupling. However, above this height while T e 
decreases monotonically with distance, T p further increases to 
2.14 x 10 6 K at 1.6 R beyond which T p declines with a slope 
steeper than T e and is overtaken by T e at 30.1 R Q . In the range 
from 1 .75 to 5 . 1 1 R Q the modeled proton temperature reproduces 
rather satisfactorily the H I measurements, indicating the heating 
mechanism we employ performs well in heating the inner solar 
wind. 

Compared with Case C, Case S yields at 1 AU a significantly 
higher proton speed and a substantially lower proton flux den- 
sity. The parameters at 1 AU now read Ve =461 km s -1 , (nv)s = 
3.14 x 10 8 cm 2 s _1 , T e = 7.5 x 10 4 K, and T p = 4.65 x 10 4 K. 
In addition, Case S yields a sonic point at rc = 3.24 R G , closer 
to the Sun than in Case C. The differences between the temper- 
atures in the two cases can be readily understood by examin- 
ing Eq.©. A higher speed in Case S leads to a more prominent 
adiabatic cooling (the second term on the left-hand-side) which 
accounts for a steeper, negative T e gradient and hence a lower 
asymptotic T e ; on the other hand, while the same also happens 
to protons, the reduced proton flux results in a higher Q^/inv) 
which leads to a steeper, positive T p gradient and eventually a 
higher asymptotic T p . However, the question remains as to how 
neglecting the curvature of the field line leads to the changes in 
the proton speed and flux? Actually this i s understandable if on e 
follows the comprehensive discussion bv lLeer & Holzerl(ll980l) : 
the more energy deposited in the subsonic region of the flow, 
the higher the proton flux; and the larger energy deposition per 
particle, the higher the terminal speed. Let us first examine the 
subsonic physics, from R to the sonic point rc- The contribu- 
tion from the electron heating, measured as an energy flux den- 
sity scaled to 1AU C c Q e (A/AE)dl, increases only slightly from 
0.3 in Case C to 0.337 in Case S (here and hereafter the energy 
flux densities are in erg cm~ 2 s _1 ). However, the effect of en- 
ergy addition to protons, which eventually derives from Alfven 
waves since L C (Q P + vF)(A/Ae)<H = F w A/Ae\^ , reduces more 
significantly from 1.089 in Case C to 0.931 in Case S. And this 
reduction is a direct consequence of the field line shape. To see 
this, we note that the wave action equation Q is a proper starting 
point, since its right-hand-side (RHS) evaluates the consequence 
of the "genuine" dissipation while its LHS involves the wave 
ponderomotive force which also contributes to the wave energy 
loss. Without dissipation, we have seen that the field line shape is 
unable to make the radial distributions different in Case C from 
those in Case S, even when the ponderomotive force is allowed 
for. Expressing Eq.© in the integral form yields 



r c a 
(1 + Ma)— — Qwaydl. 
Ae 



(7) 
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tion in the wind speed due to the field line shape is larger than 
10% throughout the latitudinal range 9 > 70°, with the largest 
relative difference being 27.8% occurring close to the equator. 
Furthermore, without further tweaking of the heating parame- 
ters, say varying the base values of wave amplitude 6vq and cor- 
relatio n length L c q from tube to tube as was done by e.g jLi et al.l 
(120041) . the modeled latitudinal dependence for 9 < 70° is 
too strong to reproduce the measured profile for the fast wind . 
Despite this, the low latitude portion of Figf3] indicates that in 
addition to a more pronounced coronal expansion rate, the field 
line curvature is also an important factor that makes the slow 
wind there slow. 



Fig. 3. Distribution with colatitude 9 of the computed solar wind 
speed v at 1 AU and the coronal expansion factor f c of the flow 
tube intersecting the Earth orbit at 9. The model result incor- 
porating (neglecting) the field line shape, called Case C (S), is 
given by the solid (dashed) line. They are obtained by varying 
the flow tube alone but keeping all other parameters unchanged. 
Cases C and S both take the radial variation of magnetic field 
strength into account, and differ only in whether the field line 
shape is considered. The blue continuous curve gives the daily 
averaged speed measurements by Ulysses during the first half of 
its latitude scan from Sep 12 1994 to Mar 4 1995. As for f c , it is 
taken to be the value of / at a heliocentric distance of 1 .7 R G . 

Now that in general dl/dr > 1, which is particularly true in the 
inner flow where Ma <k 1, the RHS of Eq.((7]) means that £) wav 
of similar amount effectively contributes more to the wave en- 
ergy reduction in the curved than in the straight case. Now let 
us discuss the asymptotic speeds. The net energy added to the 
solar wind, i.e., the electron heating plus wave dissipation plus 
enthalpy flux subtracted by radiative losses and the net electron 
heat flux flowing out of the computational domain, yields an en- 
ergy flux scaled to 1 AU of 1.495 (1.553) in Case C(S). Divided 
by the proton flux density, this gives 0.41 (0.495) xl0~ 8 erg per 
particle in Case C(S), thereby explaining the higher terminal 
speed in Case S. 

Now return to the question that, besides the tube expansion, 
what geometrical factor(s) may also affect the solar wind speed 
v. The comparison between the straight and curved cases shows 
that the field line shape could be one. In this regard, t he phy sical 
basis for the empirical relation found bv lArge et al.l d2004) may 
be simply that 9b, which measures how far the open field line 
footpoint is away from the coronal hole boundary, also charac- 
terizes the field line shape. 

What role does the shape of field lines play in shaping the 
latitudinal profile of the solar wind speed at 1 AU? This is ex- 
amined by Figure[3]which compares the speed profile in Case C 
(the solid curve) and in Case S (the dashed curve). The model 
results are obtained by varying only the flow tube along which 
the governing equations are solved, but keeping everything else 
untouched. In addition to the speed v at 1 AU, Figf3] also dis- 
plays as a function of colatitude 9 the coronal expansion factor 
f c of the flow tube reaching 1 AU at 9. Here f c is defined as 
(A/Ao)(R Q /r col -) 2 , where r col - is taken to be 1.7 R Q (see FigQJ)). 
For comparison, the solid curve in blue gives the daily aver- 
ages of Ulysses measurements of the wind speed during the first 
half of the fast latitude scan from Sep 12 1994 to Mar 4 1995. 
From Fig f3] one can see that while the speed v in both models 
is inversely correlated with f c , the results in Case C provide 
a stronger latitudinal variation. As a matter of fact, the reduc- 



4. Summary and Concluding Remarks 

Both remote-sensing measurements using the interplanetary 
scintillation (IPS) technique and in situ measurements by the 
Ulysses spacecraft show a latitudinal gradient in the wind prop- 
erties, the speed in particular, at solar minimum con dition s 
(iKojima & Kakinumdll99Cl : lMcComas et alJl200CLl2008h . What 
makes the fast wind fast and the slow wind slow still seems elu- 
sive. While there exists a robust empirical correlation between 
the coronal expansi on rate f r of the flow tu bes and the speeds 
v measured in situ (IWang&Sheelevl[T990l) . further data anal- 
ysis suggests that v depends on more than just f c dArge et al.l 
12004 . We examine whether the non-radial shape of field lines, 
which naturally accompanies any non-radial expansion, could 
be an additional geometrical factor. To this end, we solved the 
transport equations for an electron-proton solar wind, which 
incorporate the heating due to turbulent Alfven waves dissi- 
pat ed at the Kolmogorov rate, along curved field lines given by 
the lBanaszkiewicz et alJ (Il998l) model, representative of a solar 
minimum corona. The shape of field lines was found to substan- 
tially influence the solar wind parameters, reducing the terminal 
speed by as much as ~ 130 km s _1 or up to 28%, relative to 
the straight case where field line curvature is neglected. And this 
effect was interpreted in the general framewor k of energy addi- 
tion in the solar wind bv lLeer & Holzerl(ll980l) : Compared with 
the straight case, even though the wave dissipation rates may be 
similar, the field line curvature enhances the effective wave en- 
ergy deposition in the subsonic region of the flow, resulting in a 
higher proton flux and a lower terminal p roton speed . Our re sults 
suggest that the experimental finding by lArge et all (120041) may 
be interpreted in view of the fact that flow tubes having identical 
coronal expansion rates may differ substantially in their curva- 
ture. On the other hand, at solar minima even though the field 
line curvature unlikely affects the polar fast solar wind, it does 
considerably help the wind at low latitudes become slow, thereby 
helping better reproduce the Ulysses measurements. 

Could this effect have to do with the particular boundary 
conditions or the energy deposition mechanism employed here? 
This proves very unlikely. We have tested cases where the in- 
ner boundary is placed at the lower transition region where 
T = 10 s K, or the "coronal base" where T = 10 6 K. Adjusting 
the heating parameters such that realistic values can be ob- 
tained for the proton flux and asymptotic speed, and compar- 
ing solutions with field line curvature switched on and off, the 
conclusion is the same. On the other hand, using a different 
heating function where Q p is given by the dissipation of some 
unspecified mechanical energy flux at a constant scale length, 
the conclusion is the same. Moreover, this effect is not related 
to the heat conduction, even though the effective electron heat 
conductivity is reduced by a factor dr/dl in the curved case. 
Actually, turning off the electron heat flux, in which case we 
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have to fix the base density though, we reach the same con- 
clusion. All of these additional computations corroborate our 
interpretation that the changes brought forth by the field line 
shape are simply related to the differential arc-length dl involved 
in Eq.Q, whereby the more curved the field lines are in the 
subsonic region, the more wave energy reduction there, hence 
the higher the proton flux, and with the total input energy flux 
barely changing, the lower the terminal proton speed. The ef- 
fect of field line shape therefore should happen in quite general 
situations where the wave dissipation rate g wav is not directly 
proportional to the spatial derivative of background parameters 
and the energy deposition in the subsonic region plays a sig- 
nificant role in determining the solar wind flux. From this per- 
spective one may expect that the effect of field line curvature 
should not be restricted to flow tubes rooted at any particular 
latitude on the Sun. Instead, it should play a role for any flow 
tube that is curved. Moreover, this effect should also exhibit 
some influence on the flow if one employs the anisotropic tur- 
bulence treatment of the Alfven waves to incorporate the finite- 
wavelength nojiJ^XB_eflfects, such as being activ ely pursued 
at present (ICranmer et al.ll2007l IVerdini et al.l 20101) . Compared 
with a radial one, a curved line of force is likely to enhance the 
no n-WKB effect (i.e., reflectio n) in the near-Sun region (Fig. 16 
in iHeinemann & Olberd[l980l) . However, to predict what hap- 
pens next is not as straightforward as was done here. Neglect for 
the moment the wave dissipation, i.e., the waves interact with 
the flow only via their ponderomotive force. The enhanced re- 
flection tends to lead to a reduction of this force, and hence to 
a reduced mass flux as well as a lower terminal speed. Unlike 
the dissipationless case discussed in Section[2] this tendency ex- 
ists both when the wind temperat ures are prescribed (Fig. 10 in 
iMacGregor & Charbonneau| [l994) and when the win d energet- 
ics is self-consistently treated (Fig.4 in lLi & Lill2008l) . The case 
with turbulent dissipation is considerably more complicated. The 
enhancement of wave reflection in the subsonic flow means an 
enhancement of the ingoing component of the waves. With the 
turbulent dissipation trying to diminish both outgoing and re- 
flected components to similar extents, a substantial fraction of 
the reflected component may end up in the supersonic portion of 
the flow, and its consequent dissipation may lead to a higher ter- 
minal speed with mass flux barely altered. However, this is not 
the whole story since the ponderomotive force is altered in a sim- 
ilar fashion to the dissipationless case. As such the net outcome 
has to be told by a detailed numerical study, which is beyond the 
scope of this paper. Nonetheless, it seems fair to say while trying 
to reproduce the Ulysses measurements, especially those of the 
slow solar wind at low latitudes, the field line shape has to be 
accounted for. 

That said, let us say a few words on the situations where 
the curvature effects are unlikely important. This may happen, 
for instance, 1) in the spectral erosion scenario where the solar 
wind is heated by high frequency Alfven waves (with frequen- 
cies up to 10 kHz) directly launched by chromospheric magnetic 
reconnections, since now the dissipation rate is proportional to 
the directional derivative of the proton gyro- frequency and there- 
fore the background mag netic field strength (Ma rsch & Tull997t 
lHackenberg et al.l 120 00). 2) in the scenario for the solar wind 
generation where b oth the energy F w p and m ass injection F m o 



rates are fixed (e.g., Fisk 2003 ; Tu et al.ll200"5l) . To illustrate this, 
consider the scenario by Fiskl ( 20031) where the energy and mass 



injection occurs at some level where the temperature is of the 
order 1 MK, which means, roughly speaking, the contributions 
like enthalpy, heat flux and radiative losses may be neglected 
when considering the overall energy balance between the injec- 



tion point and 1 AU. Evidently the mass flux density at 1 AU is 
simply F„,qAo/Ae, and from energy conservation follows that the 
terminal proton speed is approximately ~\/2F n ,o/F m u - GM G /R Q . 
As such, the terminal mass flux density is determined by the net 
expansion of the flow tube from its footpoint to 1 AU, but the 
speed is related to neither the tube expansion nor the field line 
shape. 
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